Coherent control of trapped ions using off-resonant lasers 
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In this paper we develop a unified framework to study the coherent control of trapped ions subject 
to state-dependent forces. Taking different limits in our theory, we can reproduce two different 
designs of a two-qubit quantum gate — the pushing gate [l| and the fast gates based on laser pulses 
from Ref. |2J — , and propose a new design based on continuous laser beams. We demonstrate how 
to simulate Ising Hamiltonians in a many ions setup, and how to create highly entangled states and 
induce squeezing. Finally, in a detailed analysis we identify the physical limits of this technique and 
study the dependence of errors on the temperature. 



I. INTRODUCTION 

Trapped ions constitute one of the most promising 
systems to implement a scalable quantum computer 
In such a computer, information is stored in long-lived 
atomic states, and a universal set of gates is obtained 
by manipulating these states with lasers, and entangling 
the ions via the vibrational modes [lj. During the last 
years we have seen experimental demonstrations of var- 
ious two-qubit gates 0, H El 01 , an( i ^ remains to im- 
plement a robust scheme for scalability. Current visions 
of a scalable computer are based on the idea of moving 
the ions out of their storage area (or quantum memory) 
and make them interact in pairs, performing two-qubit 
gates U El ■ Basic steps towards the experimental imple- 
mentation of this idea have been already demonstrated. 
0. 

However, outside quantum computing there are other 
promising areas in which trapped ions may be of use. 
On the one hand, there is a great interest in preparing 
highly entangled or squeezed states, which can be used 
both for metrology |10t , or — in a more fundamental 
fashion — to characterize their decoherence properties. 
The entanglement of ions is covered in a variety of theo- 
retical papers [T^ . Il3l lldl Hlj . and it has been experimen- 
tally demonstrated for small systems 0, 0, . On the 
other hand, it has also been suggested that ion traps can 
be used to simulate a various spin Hamiltonians, ranging 
from local to long-range interactions [T^ . l2(i| . 

In all of these tasks — quantum computing, creation 
of entanglement and quantum simulation — , the goal is 
to induce some unitary evolution on the internal state 
of the ions, which is used to store the information. For 
instance, in the case of quantum computing it suffices to 
realize a phase gate between two ions 



U ph (T) 



1/4 



(1) 



while in the case of quantum simulation we want a more 



general evolution 

JJ(T) = e^'-i 5 * J *:>' ? 3+£i^*) T , 



(2) 
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where the matrices Jij and the vectors Bi determine the 
spin Hamiltonian that we want to simulate. Since the 
real interaction between the ions is described by more 
complicated Hamiltonians [See Eq. © below], any of 
these transformations is an effective one, realized after 
influencing the dynamics of the ions with external fields. 
This process, in which we dynamically modify the pa- 
rameters of the system — Rabi frequencies, detunings, 
magnetic fields, etc — , in order to achieve a well defined 
target operation, is called Coherent Quantum Control. 

While coherent control has been implicit in any pro- 
posal for quantum computing 0, H BjjlUP, El HI 
El Hi and quantum simulation pU^ofwith ion 

traps, the development of these schemes of control relied 
always on the intuition of the researcher and on a clever 
choice of approximations. In this work we show that 
many of these schemes can be translated into a unified 
framework based on state dependent forces and tunable 
traps. As characteristic examples, we will show how to 
implement the pushing gate and the fast gates based 
on laser pulses from Ref. 0. We will also propose an 
alternative and more general model of fast gate based on 
continuous off-resonant lasers. As further applications, 
we demonstrate how to induce collective H = (X^°f) 2 
and nearest neighbor H = Jcrfaf +l Ising interactions, 
and use them to produce squeezing, and generate cluster 
and GHZ states of up to 30 atoms within a extremely 
short time, T = O(l/oj), where u> is the frequency of the 
ion trap. 

Furthermore, within this unified framework we can ad- 
dress important requirements of all these coherent pro- 
cesses. Namely they should: (i) be independent of tem- 
perature (so that one does not need to cool the ions to 
their ground state after they are moved to or from their 
storage area); (ii) require no addressability (to allow the 
ions to be as close as possible during the gate so as to 
strengthen their interaction) , and (iii) be fast (in order to 
minimize the effects of decoherence during the gate, and 
to speed up the computation). All of these requirements 



2 



can be formulated as constraints of the control problem, 
and as we will see below, they can be easily solved. Last 
but not least, we study the scaling of resources as we 
try to make our control faster and answer the question 
of what is the ultimate limit for the speed of our quan- 
tum gates or entangling processes, which is shown to be 
determined both by dissipation and non-harmonic con- 
tributions to the restoring forces. 

The paper is organized as follows. In Sec. [H]we de- 
velop the formalism to study our system. First of all, 
we introduce the Hamiltonian for any number of ions 
subject to state dependent forces and quasi-lD confine- 
ment, and derive the harmonic approximation and the 
normal modes. Then, we show how to implement a uni- 
tary transformation, U = exp[i</>({erf })], made of robust 
geometrical, while leaving the motional state unchanged. 
In Sec. IIIII we apply our methods to quantum comput- 
ing in two-qubit setups. We demonstrate how to recover 
previous designs of a phase gate, including the adiabatic 
pushing gate Q and the fast gate based on n/2 laser 
pulses that kick the ions Q. Most important, since 
the generation of perfect and very short laser pulses is 
a difficult task, we design an arbitrarily fast quantum 
gate based on continuous laser sources. In Sec. IIVI we 
study the coherent control of many-ion setups. We prove 
that arbitrary spin interactions of the type Q can be 
simulated with the appropriate time-dependent forces, 
and develop a numerical method to find them. As ap- 
plications, we demonstrate numerically the creation of 
squeezed, cluster and GHZ states of up to 30 ions in a 
very short time. Sec. |V] focuses on the study of errors. 
First of all we introduce a model for dissipation on the 
vibrational degrees of freedom. Next, we solve this model 
exactly and analyze the fidelity of the effective unitary 
operations We prove that for a perfect control 

and no dissipation our scheme is insensitive to tempera- 
ture. Furthermore, the errors due to interaction with the 
environment can be computed and optimized using the 
same tools of quantum control as in Sec. lIIIlfVl Finally, 
we study the errors due to anharmonic terms in the in- 
teraction and show that both this and dissipation set an 
upper bound on the speed of the gate. In Sec. IVII we 
summarize our results and offer perspectives for future 
work. 



II. THEORY 



In this equation, V e . k is the trapping potential that con- 
fines the fc-th ion, and it may be the same for all of them 
or may change from ion to ion as in the case of micro- 
traps [l| . The time dependent external forces are denoted 
by Fi(t) and, as we will see below, they depend on the 
internal state of the ion, af . 

If we expand the previous energy around the 
equilibrium configuration without forces, given by 



(dH/d Xi ){x\ 



(0) 



,(°h - 



) = and Fi = 0, we obtain a 



set of coupled harmonic oscillators 



2m 



Fi(t) Xi 



1 

+ -jmV i jX l x.j + E . (4) 



The constant Eq is the energy of the ions at their equi- 
librium positions. The matrix V describes the restoring 
forces: it is symmetric, positive definite and it can decom- 
posed as V = MQ 2 M*, with a positive diagonal matrix 
of frequencies, £l k i = ^k$ki, and an orthogonal transfor- 
mation (MM* = M l M = I). Using this canonical trans- 
formation to define the normal modes Xi = ^ fe Mi k Q k , 



Pi 



MikPk, we arrive at the Hamiltonian 



El 

2m 



-muilQl 



, FiMikQh 



(5) 



It is now useful to write this Hamiltonian in dimensionless 
form, using the characteristic length of the oscillators, 
a k = (h/mujk) 1 / 2 , so that P k = hP k /a k and Q k = a k Q k , 
Mi k = Mi k a k . With this we arrive to 



(6) 



or, using of Fock operators, a k = (Q k + iP k )/V2, 



V2 



(J) 



B. Robust phases 

We will now demonstrate how to obtain a robust phase 
by applying forces to a harmonic oscillator. Let us begin 
with the toy model 



Normal modes 



H = fnua^a - f(t)^(^ + a) 



(8) 



As mentioned before, this paper studies a set of N 
ions, in an essentially one-dimensional confinement |39| 
and subject to some external forces. Our model Hamil- 
tonian is 



H 



E 



2m 



+ V e i(Xi) - Fi(t)Xi 



E 

i<j 



i 



47re at. 



(3) 



When integrating the Schrodinger equation associated 
with this Hamiltonian, it will convenient to use the over- 
complete basis of coherent states 



\z) :=e-^' 2 



E4=fW ='\Q + iP)- 



(9) 



Under the Hamiltonian JSJl, the coherent states behave 
somehow like classical particles in phase space, because 
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FIG. 1: Trajectories on phase space of a coherent wavepacket 
subject to a single forced harmonic oscillator, H = uaJa + 
F sin(2t)(a + oJ). In Fig. (a) we plot the usual phase space 
trajectories, (a) = {X + iP)/y/2, without forcing (dashed) 
and with F=0.1 (solid), for two initial conditions. In Fig. (b) 
we plot the same, but on a rotating frame of reference, (a) = 
{e- iut (X r +iP r ))/V2. 



their center, given by (Q) and (P), follows a classical 

trajectory, while the width of these wavepackets, given 

by the uncertainty of Q and P, remains fixed. 

More precisely, for an initial coherent state, \ip(0)) = 

\zq), the solution to the Schrodinger equation iHip = 

H(t)ip is given also by a coherent state \tp(t)) = 
p i<p(t) 



z{t)), whose phase and position satisfy 



dz 
~db 



k/(*) (* + *)■ 



dt ' 

The first equation has a solution 



z(t) = 



zo + ^fdre^m 



(10a) 
(10b) 

(11) 



that results from composing a displacement with a rota- 
tion of angular speed lo. Using the rotating phase-space 
coordinates, z r := e lut z —: Q r + iP r , we get rid of the 
motion due to the unforced harmonic oscillators and find 



2± = te^f(r)/hV2, (12a) 
d(j> dz r _ dP r dQ r dA 

The last equation means indeed that the growth of the 
phase is proportional to the area spanned by the coher- 
ent state as it moves through the phase space. The phase 
is not only geometrical in this sense, but also in the ex- 
tended definition of geometrical phase given in Ref . [2^| . 
Applying the formulas in the previous reference one finds 
that if the total phase is (f>(£), and the dynamical phase 
is always twice the geometrical one, 4>d(t) = —2(j) g (t), so 
that in the end <j)(t) — —qb g {t). 



In this work we are interested in this phase and on 
making it depend on the internal state of the particles 
governed by the oscillator. We want, however, neither 
to influence the motional state of the particles, nor to 
entangle internal and motional degrees of freedom. For 
this reason we set a time limit on the duration of the 
force and impose that after a fixed time T the coherent 
wavepacket is restored to its original state 



7(r) = 0. 



(13) 



Using this condition we derive a simple formula for the 
total phase 



4>{T) = hnj dr^- h e^f(T)z r (r) (14) 

= ^Im f 'dn / Tl dT 2 e^( Tl - T2 )/(ri)/(r 2 ). 
Jo Jo 

As a simple application, in Fig. ^ w e show the phase 
space trajectories obtained by forcing two coherent states 
with a sinusoidal force, F(t) oc sin(2uit), where oj is the 
frequency of the Harmonic oscillator. Even though look- 
ing at Fig. nja) the orbits of different initial conditions 
seem also very different, on the rotating frame of refer- 
ence the enclosed area A is always the same [Fig. db)]. 
In other words, the phase <fi is insensitive to the initial 
motional state of the system and it is thus robust. This 
property is of crucial importance when we seek applica- 
tions to real systems that cannot be cooled to the zero 
phonon limit, but which thanks to Eq. i|13|) will pick up 
the same phase regardless of the temperature. 



C. Phase of two ions 

We will now apply the results from Sec. IIIBI to a pair 
of ions. In this case there are two normal modes: the 
center of mass, x c = {x x + x 2 )/2, which oscillates with 
frequency w c , and the stretch mode, x s — x 2 — x\, which 
oscillates with frequency to s . If the ions are stored in 
the same harmonic trap, V ei k(xk) — \m<jj 2 x\ 1 these fre- 
quencies are found to be incommensurate, lu c = u) and 
uo r = cj\/3- If we store the ions in two microtraps (or 
in a more complicated potential) , the value of u r can be 
tuned from oj^/?> down to uj. 

If we exert a similar state dependent force on both 
ions, for instance by means of an off-resonance laser that 
induces a AC Stark shift on one of the internal state of 
the ions, the Hamiltonian will look as follows 

H = hjjja\a + fouj r ala r — F{t)a 7 {x\ — F(t)a 7 2Xi (15) 
= ftwa\a + huj r a\a r + F{t)(o2 — &i)d 



F{t){al+a z 2 )^={a c 



■at) 



F(t)(al-af)^=(a r + at) 



4 



where d is the equilibrium distance between the ions, 
a cr — fr/m^c r are the lengths of the oscillators, and 
of is an operator that has value +1 or —1 depending on 
whether the ion is on internal state | + 1) or | — 1). 

Since the modes are now decoupled, we can apply the 
formulas of Sec. Ill Bl almost directly. We first obtain a 
pair of commensurability relations on the force 



f dTe luJ ^ T F{r) = 0, 
Jo 



(16) 



which are just a generalization of Eq. Next we 

obtain the total phase, which up to local and global con- 
tributions is 

0=ofof/ dTl [ dT 2 G{T 1 -T 2 )F{Tl)F(T 2 ) (17) 

Jo Jo 

where Q(t) = — sin(w c |t|) — — sin(uv|i|) /2mH. 

The goal in Sec. IHII will be to tune the forces so that 
Eq. is satisfied and the phase becomes = nafa^/i, 
the required value for a two-qubit quantum phase QJ. 
We will then show two optimal solutions to this problem, 
which use either pulsed or continous forces. 



D. Phase for any number of ions 

The case of N > 2 ions exhibits a richer structure, due 
to the fact that the phase depends on all pair products 
of Oj of the polarizations of the atoms. If we apply the 
formula for the phase of one harmonic oscillator l|14|) to 
each of the modes in the effective Hamiltonian for the 
chain (0, we obtain a total phase 



i;j Jo Jo 

(18) 



with a Hermitian kernel 



M ik M jk . 

-wr sinM * l} ' (19) 



plus a generalization of Eq. (|Td|) 



V f dTe iu * T a k M ik Fi(T) = 0, Vfc. (20) 
i Jo 

In Sec. IIVI we will show how Eq. I|18|) can be related to 
an effective Ising interaction, H = J2ij jr ij cr i <T j^ whose 
precise shape can be engineered and which can produce 
interesting entangled states. 



III. FAST PHASE GATE FOR TWO IONS 

A very important application of the techniques stud- 
ied so far is to design a two-qubit quantum gate that is 



robust enough to be included in a scalable ion-trap quan- 
tum computer. This task has been pursued in a previous 
work using a slightly more complicated method, in 
which the force was achieved by kicking the ion with 7r/2 
laser pulses, and the distribution and number of these 
pulses had to be designed "by hand" . In this section we 
review this work in the light of our new formalism and 
rephrase it as an optimal control problem. This allows us 
to consider more general forces, and to find, for instance, 
a design of a gate that involves the shortest time, and 
the weakest and smoothest varying forces. 



A. Kicking forces 

We will consider two ions in a one-dimensional har- 
monic trap, interacting with a laser beam on resonance. 
The Hamiltonian modeling this system can be written as 
H = Ho + H x where Ho — fhu) c a],a c + fvjj r a\a r describes 
the normal modes of the ions and 



fi(t) 



Hi = ^ [a+e lhkxi + <r+e ihkX2 + H.c] (21) 



2 

fi(t) 



This term describes processes in which the internal state 
of an ion is changed and, as a consequence of the absorp- 
tion and emission of photons, the atom gains momentum, 
hk. The Rabi frequency f2(f) is a function of the inten- 
sity of the lasers that induce these internal transitions, 
and looking for the simplest setup we assume that it is 
the same for both ions. 

In Ref. we explained how to use the Hamiltonian Hi 
to kick the ions. The process consists on applying very 
fast laser pulses, in which the Rabi frequency, f2(i), and 

the duration of the pulse, St, satisfy J Q 0(r)dr = tt and 
St -C 27r/tt>. Let us study the evolution of the ions under 
a single laser pulse. Since the pulse is much shorter than 
a period of the trap, we can neglect the influence of Hq. 
We then use the formula 



Texp 



where II nl 



61 



, O(r) _ 
dT no 



cos(0) +ism(6)fia, (23) 



1 is a unitary vector and we impose that 

a 7r/2 pulse is produced: 9 — J Q St f2(r)/2 = tt/2. Under 
these conditions the unitary evolution is described by 



U k i 



I I iSt(liffJ+I 2 lT2) 



(24) 



If at times \t\, t 2 , . . . , t^ p } we send groups of 
{2m, 2n 2 , . . . , 2n,N p } very short laser pulses with alter- 
nating momenta, +k, —k, +k, ■ ■ ., the total operation can 
be written as U — Y[i exp[ri; x ihk{x\o\ J rx 2 <j 7 2)\, and can 
be thought of as induced by the effective force 



N p 

E 



2ni x Hka z 8{t - U). 



(25) 
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FIG. 2: (a) Trajectory in phase space of the center-of-mass 
state of the ions (Q C ,P C ) during the 2-qubit gate (solid line), 
connecting a possible initial state (black filled circle) to its 
final state (grey filled circle), for the protocol designed in 
Sec. 1111 Bl Figure (b) shows how the laser pulses (bars) dis- 
tribute in time for this scheme. Below we plot (c) the Number 
of pairs of pulses, and (b) relative angle of the two laser beams 
required to produce a phase gate using this scheme. 



performs the phase gate in a time T = 1.08/V, using 
about 4 pulses, while the second protocol allows for an 
arbitrarily short gating time T ~ N p 2 ^ 3 jv at the expense 
of using more pulses or kicks. 

The method for the first protocol is illustrated in 
Fig. [21 a-b), where we plot the phase space trajecto- 
ries followed by the center of mass mode. This se- 
quence consists basically on four groups of pulses given by 
(ni/n,U) = {(7,-ri),(l,-7a),(-l,7i ! ),(-7,n)}. The 
number n tells us how many pulses are sent within each 
group, and the parameter < 7 = cos(0) < 1.0 is 
a real number that describes how much tilted are the 
kicking lasers with respect to the axis of the trap. It 
is always possible to find a solution to Eq. i|26[) with 
ti ~ 0.538(4)(2tt/z/) > t 2 > 0. The results for the per- 
formance of the gate are summarized in Fig. ETc-d): for 
realistic values of the Lamb-Dicke parameter [5j we only 
need to apply the sequence of pulses one or two times to 
implement a phase gate. 

The second protocol performs the gate in an arbitrarily 
short time T. For shortening the time we now require six 
groups of pulses, distributed according to (ni/n,t n ) = 
{(-2, - n ), (3, -r 2 ), (-2, -r 3 ), (2, 73), (-3, r 2 ), (2, n)}, 
where the times T\, t 2 and r 3 are found numerically 
by solving the transcendental equations l|2tj|) . with the 
constraint that the whole process takes a time T = 2ri. 
As Fig. [5] shows, the number of pulses, N — 14n, 
increases with decreasing time as N p cx T -3 / 2 . This is 
just a consequence of a more general result that is shown 
later. 



When the number of pulses is odd, similar considerations 
can be done, but now a total spin flip crferf has to be 
included by hand at the end of the process. 



B. Phase gate based on kicks 

The parametrization of Eq. (|25|l is a means to simplify 
the problem of finding the optimal forces that produce 
the phase gate JU . Using the previous notation, the con- 
ditions for restoring the motional state of the ions become 



N p 



(26) 



1=1 



If these equations are satisfied, the accumulated phase 
will be 



E 

Z,m— 1 



sin(uj c ti m ) sin(uj r ti m ) 



2hk 2 nin ri 



7T/4, 

(27) 

where U m = \ti — t m \ is the time between the l-th and 
the m-th kicks. 

In Ref. we have found two possible solutions for 
these equations. The first protocol that we proposed, 



C. Phase gate based on continuous forces 

The use of 7r/2 pulses to introduce momentum in the 
ions has some inconveniences. First of all, each of the 
pulses has to be perfect, and induce a complete popu- 
lation transfer from one internal state to the other one. 
If this is not the case, systematic errors on each of the 
pulses can lead to an exponential decrease of the fidelity. 
Furthermore, as we increase the gating speed, the pulses 
may become too long to be considered as instantaneous 
kicks, and the previous formalism fails. 

What we have found, and what is also one of the main 
results of this paper, is that the phase gate may be pro- 
duced also by applying continuous forces. The search for 
this forces is then no more difficult than solving an eigen- 
value equation, where one may add restrictions such as 
smoothness of the force, and minimal total work. 

Let us take the real vector space L 2 ([0, T]) of space 
of square integrable real forces in the [0, T] interval, with 
the usual scalar product (/,<?) = J Q f(t)g(t)dt. From this 
Hilbert space, we choose a subspace TL of functions which 
are orthogonal to the modes e luJc - rt 



dr ( 



'F(r) 



dr i 



'F(r) = 0. 



(28) 
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FIG. 3: (a) Optimal forcing for a gating time Tv = 2.1 (solid), 
1.5 (dashed) and 0.9 (dash-dot), where ui = 2n x v is the fre- 
quency of the ion trap, (b) Intensity of the force vs. duration 
of the gate (solid) and visual aid (TV) -3 / 2 (dashed). All mag- 
nitudes have been adimensionalized following the text. 



Within 7i , the phase and the smoothness of the gate are 
given by cj)[F] = (F, QF) and S[F] := (F,-£sF), respec- 
tively. We will prove that the optimal (i. e. smoothest) 
force that produces a phase gate (f> — 7r/4, is simply pro- 
portional to F^, where F^ is the eigenstate 



(29) 



with largest eigenvalue \fi\. If rather than measuring the 
optimality with S[F] we use just the norm, N[F] := 
\\F\\2 — (FjF), then the eigenvalue problem is simpler 



M-^u = QF» 



(30) 



Let us prove this useful result. We have to work with 
four functionals, which are the S[F] and <fr[F] defined 
above plus two more, which measure the displacements 
originated by the force: D r>c [F] = J Q T dr e iuJ "^ T F(t). By 
choosing the space of real periodic functions which are 
orthogonal to the Fourier modes e luJa - rT we ensure that 
everything is well defined and also that the constraints 
D c ,r[F] = are satisfied. This leaves us with the problem 
of finding a force which minimizes <S[.F], while satisfying 
the last constraint <fi[F] = ir/4. There exists however a 
much easier dual problem which is formulated as finding 
the maximum of <fi[F] subject to the quadratic constraint 
S[F] = 1. Using the theory of Lagrange multipliers, this 
amounts to finding the maximum of 



C[F] = <f>[F]-n(S[F]-l), 



(31) 



where fj, is the Lagrange multiplier. Differentiating the 
Lagrangian we obtain Eq. I|29|) . from where it follows 
that /uH-Fll! — 4> is the maximal phase to be achieved, 
and the associated eigenstate F^^Jt:/{A^)/\\F\\2 is the 
force we were looking for. 

Even though we have been able to relate the control 
problem to an eigenvalue equation 129fl . there exists no 
simple analytical solution to this problem and we have to 
resort to some simple numerics. However, a very nice fea- 
ture of the two- ion problem is that, by scaling quantities 



with respect to the trap strength, w, and the wavepacket 
size, a = (h/mco)^ 1 ^ 2 , we can compute the optimal force 
independent of the setup. Using these units and expand- 
ing the force in term of Fourier modes 



/(*) = 



£ 

m=-N„ 



J2miTt/T 



F(t)a 



(32) 



we can express Eq. I|29|l as an eigenvalue equation for the 
vector c, that is to be solved numerically. The number 
of modes N m can in principle be any number above 3, 
because some degrees of freedom are lost when satisfying 
the constraint l|lrJ|) . However, the numerical experiments 
show that indeed N m — 4 provides very good solutions. 
As an example, in Fig. OJa) we show three possible forces 
for a duration of the gate Tui/2it — 2.1, 1.5 and 0.9. We 
have computed other solutions for a wider range of gate 
speeds. In Fig. 01b) we plot the mean intensity ||/|| = 
Jo 1/(^)1^ versus the total time T, and demonstrate the 
law T -3 / 2 obtained above. 

An interesting question is how much energy we have to 
put into the system in order to produce faster and faster 
gates. With our current formalism, this can be answered 
very quickly. Let us assume that we have the optimal 
force that produces a phase gate in time T <C 27r/w c r . 
Since the time is very short, we can perform a Taylor 
expansion of the function Q (r) obtaining 



Ylmh 



< 



Ylmh 



-T 3 \\F\\l 



(33) 



where ||F||i = J Q dr \F(t)\ is just a measure of the force 
applied. From here we see that 



|F||i ~ T~ 3/2 



(34) 



or, in the case of the kicked ions, N p ~ T- 3 / 2 , the scaling 
that the numerical simulations already showed. 



D. Adiabatic pushing 

As a final remark, we want to relate the methods pre- 
sented here with the pushing gate introduced in Ref . pj . 
That work proposed to trap the ions in separate mi- 
crotraps, V e ,i(xi) = ^muj 2 {xi — xf^) 2 , and to apply a 
state dependent force on two neighboring ions. This force 
should be switched on and off adiabatically with respect 
to the period of the traps, 2tt/cj, in order to approach 
the ions to each other and later on return them to their 
equilibrium positions. The adiabaticity condition ensures 
that the ions remain at all times in the ground state of the 
effective Hamiltonian which is now time dependent, 

because the equilibrium positions xf^ and the equilib- 
rium energy Eq both depend on the instantaneous value 
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of the forces. After restoring the ions to their original po- 
sitions, the only effect on the quantum state of the ions 
is a state-dependent phase 



-En 



,(°) 



(*),• 



.(0) 



(*) 



dt. 



(35) 



The previous analysis is found in Refs. P, |3(], On]. A 
very important point is that, in order not to excite the 
ions and regard the process as truly adiabatic, the forces 
have to be weak and change very slowly, and no cu- 
bic contributions to the energy should appear. In other 
words, we should be able to describe the change of E$ us- 
ing at most quadratic terms, E [x^(t)] ~ E o [x^(0)j + 

^mVijSxl dxj . Hence, rather than using the adiabatic 
theorem we can integrate the problem exactly. For a sin- 
gle harmonic oscillator we get 



*(<) 



dre^ r -/'(r), (36) 
o ^ 



where the adiabatic condition corresponds to neglecting 
the last term f'(r)/u>, and the force only has to satisfy 
f(T) = /(0) = 0. Repeating the arguments of previous 
sections, for two ions in neighboring traps the total phase 
becomes 



dr- 



2mh 



(37) 



Here lu c = lu and uj r now depend slightly on the separa- 
tion of the microtraps, but the same formula applies for 
the case in which both ions coexist in the same trap — a 
situation that could not be considered with the formalism 
of previous papers. 



IV. QUANTUM CONTROL OF SEVERAL IONS 

We will now study ID setups with more than two ions. 
As we showed before, we can still control the geometric 
phases and use them to simulate a variety of spin Hamil- 
tonians [Sec. II V A) . The design of the forces for these 
simulations is still a control problem, but a much more 
difficult one. For instance, a crucial difference is that in 
setups with more than two ions either addressability or a 
spatial modulation of the forces are required. As a pos- 
sible application of this result we study how to optimally 
generate entangled states and squeezing. In particular, 
we show that this can be done for a large number of ions 
(up to 30) in a very short time [Sec. II V B] . 



A. Simulation of spin Hamiltonians 

Given any Ising Hamiltonian 



*NxN 



and a time, T, it is possible to design a set of state depen- 
dent forces, Fi(t; J) such that after applying these forces 
for a time T , the dynamics of the ions simulates this spin 
Hamiltonian. In other words 



T 



if*H(r)dT 



-iHjT 



where T denotes time ordered product and H is the true 
Hamiltonian of the ions yp. 

The proof is very simple. Let us slice the time inter- 
val [0,T] into N 2 subintervals, < t\\ = t\ 2 < . . . < 
tjV,W-i <T. In a given time interval, / — [Uj, tij+i], we 
will switch on two forces, and leave all other ions on their 
own, 

F k (t) = 0, t € [Uj, ti,j+i], Vfc ^ 

The active forces Fi(t) and Fj(t) must satisfy several 
equations 



= J e WkT M ka F a {T)dT, a = i,j, k=l...N (38) 



dridr 2 Fi(ri)Gij(n - T 2 )Fj(T2). 

i J i 

It is not difficult to convince oneself that these equations 
always have a solution, and that by repeating this proce- 
dure we will get an effective total phase, </>, that resembles 
the one produced by the Ising model during a time T. □ 

We have to make several remarks here. The first one 
is that since the operator that we want to simulate is 
symmetric, Jy = Jji, and since the diagonal terms only 
contribute to a global phase, the number of intervals can 
be actually decreased to N(N — l)/2. 

However, more important is the fact that we can use 
coherent control to find optimal forces, F, which instead 
of piecewise continuous are the smoothest possible and 
have the optimal norm, while giving rise to the same 
effective Hamiltonian. This task has been performed nu- 
merically for some models, and the results will be shown 
in the following section. 

From the point of view of quantum simulation, we 
would like to be able to model more than just an Ising 
model, which is essentially classical. For instance, one 
would like to be able to introduce transverse magnetic 
fields, Yli a i°~f> 01 to simulate a Heisenberg interaction, 
(JiGj, and in general, a unitary operation of the form iffil 
would be desirable. The answer to this problem is once 
more the stroboscopic evolution, or a Trotter expansion 
of the operator @, 



U 



n 




(39) 



In this expansion, we decompose the total unitary as a 
product of phase gates, that are originated by forces that 
depend on of, of and of . In practice, one would switch 
on a magnetic field h\ and perform a phase gate with 
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coefficients Jf- for a time T/N, rotate the spins so that 
Uy becomes a z , apply the phases with Jfj, etc. 

It is also worth noticing that if we switch on the state- 
dependent forces acting on different ions, and make them 
oscillate with a single frequency Q around a constant 
value, Fi(t) = /i sin(Slf)of , for a long time, the effective 
interaction is a particular spin Hamiltonian 



ijk 



M ik M jk [l 



Amh(ujl - fi 2 ) 



' )li ° ] fiO*o* + oA), (40) 



In the limit — > 0, this model corresponds to the one 
found in |l9j. As it was shown there, depending on 
whether the forces operate longitudinally or transversely 
to the ion trap, this continuous force will give rise to long 
range or short range interactions. 




FIG. 4: Generation of a (a,c) cluster state of TV = 10 ions 
and of (b,d) a GHZ state of N = 20 ions, using common 
forcing for a time T = l.l/u>. (a,c) Time dependence of the 
forces, /(f). (b,d) Fidelity with respect to the target state as 
a function of time. 



B. Coherent control and design of entanglers 

The simulation of an Ising interaction is by itself inter- 
esting, and has important applications such as creation of 
many-qubit quantum gates, quantum simulation and adi- 
abatic quantum computing. However, a most straight- 
forward and useful application of an Ising Hamiltonian 
is the generation of many-body highly entangled states 
called graph states 32]. Roughly speaking, let us imag- 
ine that we have a set of N spins, which we represent 
by points or vertices, and let us connect these points by 
lines or edges. The resulting graph can be described by 
an adjacency matrix which takes value = 1 only if the 
spins i and j are connected. To each graph we can thus 
associate a Hamiltonian of the form Hj = J2i j Jij^i^j- 
It has been shown that after applying this interaction 
over a certain time on a transversely polarized state, the 
outcome is a highly entangled state called a graph state: 

\TpG) = ^e- i % H °(\0) + \l)f N . (41) 

When the graph has a latti ce g eometry, these states are 
also known as cluster states |33j , and form the basic ingre- 
dient of the one-way quantum computer |.3 l|. However, 
a particularly important case without lattice geometry is 
the GHZ state, 

|GHZ> ~ |0>® iV + |l>® JV , (42) 

which is essentially generated by the interaction Jy = 1 
or Hj = (J2i of) 2 . The GHZ state is one of the best stud- 
ied entangled states, it constitutes a canonical example of 
Schrodinger cat state, and it could have important appli- 
cations in the field of precision frequency measurements, 
providing a 1 /yN precision increase for N entangled ions 
|l0l ITU, a point already demonstrated experimentally in 

Ref. nEa. 

We have investigated how to implement these highly 
entangled states using our quantum control techniques. 



The idea is very simple: we design a matrix for 
our graph state, and look for the time- and state- 
dependent forces that implement the phase transforma- 
tion exp(— inHjir/A) within a fixed time T. For sim- 
plicity, even though it is not warranted to succeed, we 
look for forces that have a common time dependence 
Fi(t) = Jif(t)af, 1 7(| < 1. These forces could be imple- 
mented with an appropriate intensity mask, which deter- 
mines the relative amplitudes ji, and a global intensity 
modulation, which gives the function /(f). Expanding 
this modulation in Fourier modes, we find N(N—l)/2+N 
equations which define a possible entangling procedure. 
We have solved numerically these equations, both for the 
GHZ state and for the cluster state. While in the first 
case we always found exact solutions with a small num- 
ber of modes (i.e. 50 modes for 30 ions), the generation 
of the cluster state was always approximate with high 
fidelity, F ~ 99.9%. The error in this case has its origin 
in our particular choice of forces. 

In Fig. 0] we show the entangling procedure for a setup 
with 10 and 20 ions, even though chains of up to 30 ions 
have been considered. We measure the fidelity of the pro- 
cess as the overlap between the target state and the one 
achieved. If 6 J is the difference between the desired inter- 
action and the achieved one, then F = ^tv Yls e ~ lS SJTs , 
where the sum is performed over all possible spin con- 
figurations, Sk = ±1. The time scales for the generation 
of the interaction are independent of the size of the sys- 
tem, and for instance we can produce a GHZ state of 
20 ions in a time T = 1.1/w, to be compared with the 
time T = 30007r/a-> required when individually address- 
ing one of the vibrational modes 0]. The strength of 
the forces, though, grows moderately with the number of 
ions, which can be inconvenient. However, thanks to the 
periodicity of the forcing, /(f), if we divide the intensity 
of the forces by a factor of 2, /(f) — > /(f)/2, the same 
gate be produced in a longer time AT. Furthermore, the 
forces that we present in this paper are not optimal, and 
have been found with a straightforward Gauss-Newton 
method. If high fidelity is not required, one may find 
better solutions with fewer modes, but most important 
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we expect significant improvements by the application 
of better numerical algorithms to search for the optimal 
entanglers. 

Using the Ising interaction we can not only produce 
graph states, but also squeezed states: states in which the 
variance of one spin component, AS X , has been decreased 
at the expense of increasing the other variances. As it 
was shown in l35| . a Hamiltonian of the form H = 
(single axis squeezing) or H = J% — Jy (two-axis squeez- 
ing) can be used to produced squeezing. Both models 
can be simulated using our tools, either directly, as in 
the single-axis squeezing, or stroboscopically, for the XY 
interaction. Indeed, the stroboscopic simulation of the 
two-axis squeezing resembles the scheme of tt/2 pulses 
used in [36l | to effectively switch off the interaction in 
two-mode Bose-Einstein condensates |37j . 

V. OPTIMAL CONTROL OF ERRORS 

Up to now, we have assumed that the motion of the 
ions is not disturbed during the time when the controlling 
forces are applied. In this section we will show how to 
take these effects into account for a realistic model of 
dissipation. The main result is that the fidelity of the 
process can still be computed and that the there are two 
sources of error: one due an imperfect control of the ions, 
which introduces some temperature dependence on the 
fidelity [Sec. IV (J] , and another one due to the dissipation, 
that can be treated as another constraint for the control 
problem [Sec. IV B] . Finally, we will comment on possible 
extensions outside the harmonic regime [Sec. IV D| 

A. The model and a exact solution 

We study the dissipation with a master equation that 
arises from coupling the phonon modes with a "classical" 
bosonic bath in thermal equilibrium 

~d± P = ~h^ Hl,P ^ + ^2^ N k( 2a {p a k - ako[p - pa k a\) 

+ ^2lk{N k + l){2a k pa\- a\a k p- pa\a k ). (43) 
k 

Here j k describes the coupling to an external bath and 
N k is the mean number of bosons on that bath and it is 
related to its temperature. The Hamiltonian in Eq. 143(1 
is written in the interaction picture 

Hi = + <7fci4K. (44) 

ki 

In order not to obscure the discussion, we will assume 
that each phonon mode interacts with an independent 
bath. In that case, the forces in the rotating frame of 
reference become, hg k i{t) — -^Fi(t)M ik a k e~ tuJkt . How- 
ever, it is easy to generalize the following analysis to a 



more realistic model in which each ion couples indepen- 
dently to the environment, and the operators a k and a k 
do not represent the phonons, but the displacements of 
the ions. 

To study the fidelity of a gate we only need to know 
the matrix elements of the reduced density matrix for the 
internal degrees of freedom. This matrix may be written 
as a collection of expectation values, 

Preal = ^ (E s ,r) £ S ,r , (45) 
s.r 

where S SjI . = \si . . . s N )(ri . . .r N \, r k ,s k = ±1 forma 
complete basis for the space of 2 N x 2 N complex matrices. 

The calculations that provide us with the evolution of 
(E s , r ) are detailed in Appendix ^ Here we will only 
summarize the main result, which is that the reduced 
density matrix can be written as 

(E. T )(t) = e -"->*e-*£j* J ^ s *(£ s P )(0)e*£i* J i* r * p *. 

(46) 

In other words, the spin density matrix has the form 

Preal(t)=e- C [U Preal {0)U^], (47) 

where U = exp(— Jjk&j& k ) is the operation that we 
want to perform, and L(p) = J2 S r K s,rE Sir /?£ SiI . is re- 
sponsible for the decay of coherences. 

In comparison with the previous part of the paper, the 
unperturbed orbits in phase space, that is the evolution 
without external forcing, are now not circular orbits, but 
by spirally decaying ones. This fact translates in the new 
conditions for uncoupling internal and motional degrees 
of freedom lj2U|) 

/ dTe~ [iu > k+lk)T F,j{t) = 0, Mj,k (48) 
Jo 

which now depend on the exponential decay rate, 7^, 
given by our dissipation model. This model-dependence 
is also evident in the kernel that produces our unitary 
operation, Jy , which now reads 

M^E^^MV-'''. (49) 

k 

Finally, we would like to remark that the conditions, the 
phase and the kernel are only slightly modified when we 
consider a local coupling to the environment. 

B. Quantum control of errors due to dissipation 

To understand better the implications of Eq. (|47jl , let 
us study the dynamics of a single ion. In this case the re- 
duced density matrix can be expressed uniquely in terms 
of the expectation values (a z ) and (c + ). Furthermore, 
since the magnetization is constant, we can compute the 
Uhlmann fidelity exactly as a function of (u + ). Let us 
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assume that initially the system is in a pure state, and 
define (a + } rea i — e lS ^ K (a + ) i d, where the subindex "id" 
denotes the ideal value obtained in absence of errors. The 
Uhlmann fidelity of the gate is [38| 



HpreaUPid) = a/1 + |(<r+)id| 2 (2e-«cos(<^) - 1). (50) 

Two types of errors contribute to the decay of the fi- 
delity. The first type is made of control errors. These 
errors contribute both to the spurious phases (6(j>) and 
to an effective decay due to not restoring the motional 
state of the ions (k =t because (3{T) ^ and U(T) ^ 1). 
These errors cause a smooth dependence on the temper- 
ature to appear, as shown later in Sec. IV CI 

The second type of errors are due to dissipation dur- 
ing the forcing of the ions. Their contribution to the 
exponential decay is 



ISSIp 



= 7(AT+|) / rfr|/3(r)p 
'o 



(51) 



One would be led to think that if dissipation acts on a 
much larger scale than the time required to perform our 
gate we can neglect it completely. However, a simple ar- 
gument shows that this is not the case. As we saw before 
in Sec. IIII CI the strength of the forces scales roughly as 
F ~ T~ 5 . This scaling allows us to give a worst case 
estimate of (3(t) ~ Ft and conclude 



issip 



1 T- 



T < 1/w. 



(52) 



What this means is that slower gates will involve smaller 
displacements of the ions, which in turn translates into 
less dissipation. On the other hand, a too long applica- 
tion of a force also gives more time for the dissipation 
to operate and the optimal duration should be a com- 
promise between both processes. It is thus possible and 
recommended to optimize the forces F(t) taking not only 
into account the properties of the force (i.e. differen- 
tiablity and intensity), but also trying to minimize the 
decay k induced by the force. From the numerical point 
of view, this new control problem is only slightly more 
complicated than the ones we have solved in Sees. IIIIIIVI 
because Kdissip is a nonlinear function of the forces. 



plus a perturbation of the phase 



(54) 



which can be interpreted as a change in the effective in- 
teraction between ions [Sec. HYj 



SJij = dndT2Qij{T2 — n) x 

Jo Jo 

We will assume as initial condition a pure state of the 
internal degrees of freedom and a thermal state of the 
vibrational ones p(0) = \tp}(ip\ <S> p V ib{T). The Uhlmann 
fidelity at the end of the process is given by 

F{Pid,Preal{t)) = Tr vib ((^\U^ d p(t)U ld \lp)) 

= (ip\Tr vib [ulp(t)U id ]\i)). (55) 
Expanding \tp) — c s \si . . . sn), we obtain 

s.r 

When we neglect dissipation, the previous expectation 
values can be computed in terms of the final displace- 
ments, (3k(t), and the residual phases as follows: 



'Uj dP (t)U, 



d =Tr|n^[/3 fe (t)( Sfe -r fe )]p OTb (r)| ^57) 



xe 



(s'Hs-r'Wr) * 



Here D(z) = exp(za^ — z*a) is the displacement operator, 
Pvib{T) is a thermal state 



JV 



P mb {T) = ($^p- 1 1 dzdze-^*/ k ° T \z) k (z\, 



k=l 

(58) 

and thus (D(z)) = exp[-\z\ 2 (l/2+k B T/huj k )} = exp(C fe ) 
so that the total fidelity becomes 

T = J2 \c s \ 2 \\c r \ 2 e listSJs ~ rt5Jr) e Ck(s ^ 2 . (59) 



C. Errors due to an imperfect control: influence of 
temperature 

Let us denote by Uid — exp(^- J Jijcrfaj) the ideal 
operation that we want to perform, and by U rea i the op- 
eration with errors. In this section, the only source of 
error that we consider is an imperfect control, denoted 
by a perturbation of the state-dependent force induced 
on the ions, (Fi(t) + 6Fi(t))a^ . According to the previous 
analysis, the effect of this perturbation will be a residual 
state-dependent displacement of the coherent wavepack- 
ets at the end of the process, f3k(T), 



0h = ~i 



dr i 



kT a k M ik 5Fi{t) 
1 = 1 



=:^/3 fci a?. (53) 



D. Errors due to larger displacements 

The previous studies can be generalized to arbitrary 
interactions and trapping potentials. Let us assume a 
complicated Hamiltonian 



+ V(xi,...,x N ), (60) 



describing the traps and the ion-ion interaction. The 
evolution equation for the position of the ions are of the 
form 



k,i 



. _pk . 

%i 7 Pi 

m 



dV 



fi(t)o-l 



(61) 
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Since the operators of are conserved quantities, the pre- 
vious equations can be thought of as a simple problem 
of Newtonian mechanics, even though in practice, both 
Xi and pi are operators. We can thus represent a general 
solution as {xi(t; a),pi(t; a)}, where a denotes the val- 
ues of of operators. The phase of the ions is then com- 
puted by analyzing the evolution of the operators. 
These operators must undergo a unitary transformation 
<J + {t) — U (t)pf (0)W (t) in which the dependence on the 
{erf} operators must be of the form |4[j 



U(t) = exp[]T Biifiot + J2 <M*K ? 



(62) 



F max T 2 /2m, where F max is the maximal force applied on 



the ions. Next we will use the scaling i 



2 T 5 F max /hm 



to show that roughly (x max /a) ~ 4>/AujT. With this, and 
first order perturbation theory we compute the error and 
estimate it as 



Fanh — 



(§)" 



A3/2 



4 3 / 2 u;T' 



(68) 



If we want to apply our phase operations to build a quan- 
tum computer, we need E < 10~ 4 and there is a limit 
on the speed of the gate T > 10~ 3 /a->, which nevertheless 
gives gating rates of the order of 100 MHz. 



Using the commutation relation [af ,<7j] = —25ijafcrf, 
we obtain 

a+(t) = e Ml *f +2 £j**«(*K*J .+ (o). (63) 
Combining this with the Heisenberg equation for af 

^=2 M *)^=2^^ ) (64) 

we find, up to global phases, 

i-T 



I drfi{j)xi{T)al 
i Jo 



(65) 



From this analysis we see that we must impose two 
conditions on the process. On one hand, the orbits of 
the ions must be periodic so as to disentangle the internal 
and motional degrees of freedom 



Xi(T;a) ~ Xi(0), Pi(T;a) ~ Pi(0). 



(66) 



On the other hand, the phase <j) must be independent 
of the initial conditions, {xj(0),pj(0)}. Satisfying both 
conditions is impossible in general, but if we restrict 
ourselves to small displacements and harmonic restor- 
ing forces, dV/dxi ~ J^j Vij x j> it is possible to integrate 
Eq. H61fl and recover our expressions for the phases (|18H . 

If, however, the qubit and higher terms in 
V(xi, . . . , xn) become important, we will fail the restor- 
ing condition l|66|l . and induce some entanglement be- 
tween the motion and the spin of the ions. The errors 
due to these anharmonic terms are of the order 



d 3 V 



dxidxjdxi 



XiXj X{ 



3Hu>a 



(67) 



where a ~ ^Jh/raoj is the length associated to the har- 
monic oscillator, d 3 — e 2 ^irevmuj 2 is the equilibrium 
distance between two ions, and a; is a typical displace- 
ment. Since a trivial analysis of these errors is not possi- 
ble, we can only produce a pessimistic, first order bound 
that restricts the error induced by this perturbation on 
the wavefunction. First we will give a worst case predic- 
tion for the maximum displacement of the ions as x max < 



VI. CONCLUSIONS 

We have developed a unified framework to study the 
coherent control of trapped ions by means of state- 
dependent forces and robust geometric phases. Our tech- 
niques can be used to perform fast two-qubit gates be- 
tween pairs of ions. For an adiabatic switching of the 
forces and for the case of pulsed lasers we are able to 
reproduce the proposals of and Q, and with very lit- 
tle work we can design the optimal forces that produce 
a phase gate in a given time with the lowest intensity. 
Using the same tools and a larger number of ions, we can 
simulate either continuously or stroboscopically a num- 
ber of spin s = 1/2 Hamiltonians. Furthermore, we are 
also able to create highly entangled states and squeez- 
ing, and as prototypical examples we have shown how 
to obtain a GHZ state of 20 ions in a very short time, 
T = 1.1/uj. Finally we have studied the sources of error 
in the application of our gate, which are an imperfect 
control, dissipation and anharmonic terms in the inter- 
action. The first type of errors could be ideally corrected 
and introduce a smooth decay of the fidelity with the 
temperature. The second type of errors induces also a 
decay of the fidelity, but the amount of this error can 
be optimized using the tools of quantum control. Both 
dissipation and anharmonicity set up upper limits on the 
speed of the gate. This limit is however very weak, since 
it allows theoretically a gating speed of hundreds of MHz, 
and it could be overcome by a numerical study of the role 
of anharmonic terms in the motion of the ions. 

While concluding this paper we became aware of the 
work by P. Staanum, M. Drewsen and K. M0lmer [2^ on 
performing quantum gates using continuous laser beams. 
The ideas shown in Ref. [2^ are equivalent to the de- 
velopment of a two-qubit gate done in Sect. , with the 
difference that we provide an optimal solution for the 
problem. 
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APPENDIX A: SOLUTION OF THE MASTER 
EQUATION 

As we mentioned before, the density matrix is charac- 
terized by the expectation values (T, s r ). However, it is 
much more difficult to work with these expectation val- 
ues, than with 



.4 



E. )P V(t) := S s 



(Al) 



By imposing that V(0) = 1 and that V(T) is at most a 
phase, we will be able to relate the reduced density ma- 
trices preai(0) and Preai(T). It is easy to see that indeed 
the operator V(t) is a displacement operator and that 
Eq. (|A1|) is essentially the solution of the nondissipative 
case, where (3 k measures the separation in phase space 
between configurations with internal states s and r. 

The equation for the expectation value of an arbitrary 
operator, A, can be written as follows 



dt 



(A) = {a t A) + j i {[H I ,A])-'y((D a A)a + at{D a )) 



r yN(D a D a i-A + D a iD a A). 



(A2) 



Here, D a A ;= \A,<v] and D a \A := [a, A] are two su- 
peroperators which first of all commute, [D a ,D a t] = 0, 
and second they are related to the formal derivatives 
with respect to the operators a and a) . So, for instance, 
D a f(a,a ji ) = d a f(a,a i ) for any analytical function /. 
If we substitute Eq. IjAlfl into Eq. (|A2|) , and use 

„ . v-^ 9 A ■ OA ■ . 



= u 



[M - §/3£) +fa(-a k ~ , (A3) 



we will obtain 

k 

+ E ( A [M - + km - M) 

k 

+ J2^N k (A(-l3ta k + (3kai)) 

k 

- ^ 7 /c(^ + i)|/3 fe | 2 <A), (A4) 



whith the new parameters g kr '■= J2i 9ki{t)ri. 

Here is where we impose a particular evolution of the 
displacements on phase space, $k+lk/3k+i(gk r ~ 9ks) = 0- 
This equation has a trivial solution 

Pk{t)=i f dre-^-^[gi T {r)-gt s {r)], (A5) 

After substituting this value all terms containing Fock 
operators are cancelled and we are left with 



dt 



(A) 



-k + i<j>)(A), 



where the decay K is 



K(*) = y)7fc(^fc + |) / dr|&(r) 

u JO 



and the total phase <ft — Jij( r i + s i)( r j 
tcrmined by the matrix 



(A6) 

(A7) 
sj) is de- 



ft /-Tl 

Jij—Imj dn I dT2gki{n)g k j{T<z)*e 







* p -7*>(t'2-ti) 







(A8) 

Using the symmetry of this matrix, the formula for the 
phase can be rewritten as <f> = J2a Jij{ r i r j ~ s i s j), an d 
the results mentioned in Sect. IV Al quickly follow. 
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This is only for illustrative purposes. More complex se- 
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We use the identity (erf) 2 = 1 



